This file reads the log likelihood of all pairs apart same distance of the simulated datasets.

Now plot the variance of lnL with varying log_p in each distance

for (sim in 1:100){
  grid.lnL <- get(paste("Tract_", tract, ".0_sim_", sim, "_grid_lnL", sep = ""))
  sim.info <- get(paste("sim.tract.", tract, sep = ""))
  lnL.var <- apply(grid.lnL, 2, var)
  plot(lnL.var, type = "l", 
     xlab = "distance", ylab = "variance in lnL",
       main = sim)
  plot(exp(-grid.lnL[, 1]), rowSums(grid.lnL[, -1]), type = "b",
       xlab = "Tract length", ylab = "lnL", 
       main = paste("lnL vs Tract sim = ", sim))
  print(paste("In simulation, number of IGC events:", sim.info["num IGC", paste("sim", sim, sep = "_")], 
              ", mean potential tract length:", sim.info["mean potential tract length", paste("sim", sim, sep = "_")],
              ", mean realized tract length:", sim.info["mean realized tract length", paste("sim", sim, sep = "_")]))
  print(paste("Max lnL at grid value:", round(exp(-grid.lnL[, 1])[which.max(rowSums(grid.lnL[, -1]))])))
  if (paste("sim", sim, sep = "_") %in% colnames(PSJS_HKY_Tract_50_combined_summary)){
    print(paste("PSJS estimate:", PSJS_HKY_Tract_50_combined_summary["tract_length", paste("sim", sim, sep = "_")]))
    
  }
}

## [1] "In simulation, number of IGC events: 22 , mean potential tract length: 53.0909090909091 , mean realized tract length: 53.0909090909091"
## [1] "Max lnL at grid value: 15"
## [1] "PSJS estimate: 17.2264286677101"

## [1] "In simulation, number of IGC events: 29 , mean potential tract length: 70.4137931034483 , mean realized tract length: 58.5172413793103"
## [1] "Max lnL at grid value: 25"
## [1] "PSJS estimate: 25.8290905925275"

## [1] "In simulation, number of IGC events: 25 , mean potential tract length: 45.4 , mean realized tract length: 39.52"
## [1] "Max lnL at grid value: 55"
## [1] "PSJS estimate: 55.0163671864766"

## [1] "In simulation, number of IGC events: 31 , mean potential tract length: 57 , mean realized tract length: 54"
## [1] "Max lnL at grid value: 40"
## [1] "PSJS estimate: 40.0806223939151"

## [1] "In simulation, number of IGC events: 27 , mean potential tract length: 53.3703703703704 , mean realized tract length: 52.037037037037"
## [1] "Max lnL at grid value: 70"
## [1] "PSJS estimate: 66.9784593008133"

## [1] "In simulation, number of IGC events: 28 , mean potential tract length: 63.1428571428571 , mean realized tract length: 55.1785714285714"
## [1] "Max lnL at grid value: 100"
## [1] "PSJS estimate: 115.93519145024"

## [1] "In simulation, number of IGC events: 33 , mean potential tract length: 44.0606060606061 , mean realized tract length: 38.5454545454545"
## [1] "Max lnL at grid value: 30"
## [1] "PSJS estimate: 29.3227195815801"

## [1] "In simulation, number of IGC events: 27 , mean potential tract length: 61.7407407407407 , mean realized tract length: 60.962962962963"
## [1] "Max lnL at grid value: 35"
## [1] "PSJS estimate: 34.9817003225787"

## [1] "In simulation, number of IGC events: 31 , mean potential tract length: 49.258064516129 , mean realized tract length: 46.6129032258064"
## [1] "Max lnL at grid value: 25"
## [1] "PSJS estimate: 27.0089860007074"

## [1] "In simulation, number of IGC events: 30 , mean potential tract length: 40.1 , mean realized tract length: 39.5666666666667"
## [1] "Max lnL at grid value: 55"
## [1] "PSJS estimate: 54.6869850803605"

## [1] "In simulation, number of IGC events: 26 , mean potential tract length: 48.7692307692308 , mean realized tract length: 31.8461538461538"
## [1] "Max lnL at grid value: 20"
## [1] "PSJS estimate: 18.2838640824187"

## [1] "In simulation, number of IGC events: 31 , mean potential tract length: 37.1612903225806 , mean realized tract length: 35.7741935483871"
## [1] "Max lnL at grid value: 15"

## [1] "In simulation, number of IGC events: 27 , mean potential tract length: 41.8888888888889 , mean realized tract length: 39.3703703703704"
## [1] "Max lnL at grid value: 35"
## [1] "PSJS estimate: 32.8752954509935"

## [1] "In simulation, number of IGC events: 24 , mean potential tract length: 46.5 , mean realized tract length: 39.6666666666667"
## [1] "Max lnL at grid value: 30"
## [1] "PSJS estimate: 28.2827425528188"

## [1] "In simulation, number of IGC events: 29 , mean potential tract length: 40.9655172413793 , mean realized tract length: 40.5862068965517"
## [1] "Max lnL at grid value: 40"
## [1] "PSJS estimate: 39.8980249835229"

## [1] "In simulation, number of IGC events: 20 , mean potential tract length: 49.35 , mean realized tract length: 45.6"
## [1] "Max lnL at grid value: 40"
## [1] "PSJS estimate: 40.2967432873549"

## [1] "In simulation, number of IGC events: 38 , mean potential tract length: 41.2631578947368 , mean realized tract length: 35.4473684210526"
## [1] "Max lnL at grid value: 45"
## [1] "PSJS estimate: 45.8862886128301"

## [1] "In simulation, number of IGC events: 34 , mean potential tract length: 55.7941176470588 , mean realized tract length: 52.8529411764706"
## [1] "Max lnL at grid value: 35"
## [1] "PSJS estimate: 34.0672683704242"

## [1] "In simulation, number of IGC events: 23 , mean potential tract length: 60.4347826086956 , mean realized tract length: 58.2608695652174"
## [1] "Max lnL at grid value: 20"
## [1] "PSJS estimate: 17.6005712156941"

## [1] "In simulation, number of IGC events: 23 , mean potential tract length: 57.304347826087 , mean realized tract length: 50.9565217391304"
## [1] "Max lnL at grid value: 35"
## [1] "PSJS estimate: 35.1509619137355"

## [1] "In simulation, number of IGC events: 28 , mean potential tract length: 55.1071428571429 , mean realized tract length: 42.9642857142857"
## [1] "Max lnL at grid value: 10"
## [1] "PSJS estimate: 12.2006430230297"

## [1] "In simulation, number of IGC events: 33 , mean potential tract length: 53.0606060606061 , mean realized tract length: 44.1212121212121"
## [1] "Max lnL at grid value: 35"
## [1] "PSJS estimate: 34.3586355344462"

## [1] "In simulation, number of IGC events: 27 , mean potential tract length: 58.8148148148148 , mean realized tract length: 49.5555555555556"
## [1] "Max lnL at grid value: 70"
## [1] "PSJS estimate: 65.3126299196728"

## [1] "In simulation, number of IGC events: 33 , mean potential tract length: 44.3030303030303 , mean realized tract length: 40.8181818181818"
## [1] "Max lnL at grid value: 25"
## [1] "PSJS estimate: 24.6704111417318"

## [1] "In simulation, number of IGC events: 26 , mean potential tract length: 64.8461538461538 , mean realized tract length: 57.1538461538462"
## [1] "Max lnL at grid value: 35"
## [1] "PSJS estimate: 37.1173902200211"

## [1] "In simulation, number of IGC events: 28 , mean potential tract length: 58.6071428571429 , mean realized tract length: 55.8928571428571"
## [1] "Max lnL at grid value: 40"
## [1] "PSJS estimate: 38.3774116684469"

## [1] "In simulation, number of IGC events: 30 , mean potential tract length: 62.8 , mean realized tract length: 46.3"
## [1] "Max lnL at grid value: 30"
## [1] "PSJS estimate: 31.5257597396145"

## [1] "In simulation, number of IGC events: 23 , mean potential tract length: 48.0434782608696 , mean realized tract length: 41.7391304347826"
## [1] "Max lnL at grid value: 30"
## [1] "PSJS estimate: 32.1741563324853"

## [1] "In simulation, number of IGC events: 16 , mean potential tract length: 42.3125 , mean realized tract length: 42.3125"
## [1] "Max lnL at grid value: 15"
## [1] "PSJS estimate: 16.4247346499549"

## [1] "In simulation, number of IGC events: 26 , mean potential tract length: 55.3846153846154 , mean realized tract length: 53.6538461538462"
## [1] "Max lnL at grid value: 45"
## [1] "PSJS estimate: 46.9870580229305"

## [1] "In simulation, number of IGC events: 32 , mean potential tract length: 68.53125 , mean realized tract length: 61.46875"
## [1] "Max lnL at grid value: 40"
## [1] "PSJS estimate: 39.0402804830966"

## [1] "In simulation, number of IGC events: 32 , mean potential tract length: 54.5 , mean realized tract length: 49.5625"
## [1] "Max lnL at grid value: 30"
## [1] "PSJS estimate: 31.038360322701"

## [1] "In simulation, number of IGC events: 39 , mean potential tract length: 55.5128205128205 , mean realized tract length: 53.8461538461538"
## [1] "Max lnL at grid value: 35"
## [1] "PSJS estimate: 34.0788875123054"

## [1] "In simulation, number of IGC events: 22 , mean potential tract length: 44.2727272727273 , mean realized tract length: 34.5"
## [1] "Max lnL at grid value: 15"
## [1] "PSJS estimate: 17.2055374114319"

## [1] "In simulation, number of IGC events: 26 , mean potential tract length: 51.6538461538462 , mean realized tract length: 45.5769230769231"
## [1] "Max lnL at grid value: 25"
## [1] "PSJS estimate: 26.5111712118757"

## [1] "In simulation, number of IGC events: 24 , mean potential tract length: 62.25 , mean realized tract length: 36.625"
## [1] "Max lnL at grid value: 50"
## [1] "PSJS estimate: 48.3867529956448"

## [1] "In simulation, number of IGC events: 27 , mean potential tract length: 53.4074074074074 , mean realized tract length: 47.5925925925926"
## [1] "Max lnL at grid value: 15"
## [1] "PSJS estimate: 17.1216108633329"

## [1] "In simulation, number of IGC events: 29 , mean potential tract length: 40.7931034482759 , mean realized tract length: 34.6896551724138"
## [1] "Max lnL at grid value: 20"
## [1] "PSJS estimate: 21.2594791663957"

## [1] "In simulation, number of IGC events: 23 , mean potential tract length: 44.1304347826087 , mean realized tract length: 42.9130434782609"
## [1] "Max lnL at grid value: 25"
## [1] "PSJS estimate: 23.5819112551754"

## [1] "In simulation, number of IGC events: 29 , mean potential tract length: 45.9310344827586 , mean realized tract length: 40.5172413793103"
## [1] "Max lnL at grid value: 40"
## [1] "PSJS estimate: 38.2799724162543"

## [1] "In simulation, number of IGC events: 32 , mean potential tract length: 69 , mean realized tract length: 59.78125"
## [1] "Max lnL at grid value: 50"
## [1] "PSJS estimate: 49.2399603523822"

## [1] "In simulation, number of IGC events: 35 , mean potential tract length: 56.2857142857143 , mean realized tract length: 44.1714285714286"
## [1] "Max lnL at grid value: 50"
## [1] "PSJS estimate: 50.715379806399"

## [1] "In simulation, number of IGC events: 21 , mean potential tract length: 50.7619047619048 , mean realized tract length: 48.3333333333333"
## [1] "Max lnL at grid value: 25"
## [1] "PSJS estimate: 25.526137294565"

## [1] "In simulation, number of IGC events: 27 , mean potential tract length: 42.5185185185185 , mean realized tract length: 40.6666666666667"
## [1] "Max lnL at grid value: 25"
## [1] "PSJS estimate: 22.865296492081"

## [1] "In simulation, number of IGC events: 27 , mean potential tract length: 54.5185185185185 , mean realized tract length: 52.7037037037037"
## [1] "Max lnL at grid value: 25"
## [1] "PSJS estimate: 24.4462162891691"

## [1] "In simulation, number of IGC events: 22 , mean potential tract length: 40.8181818181818 , mean realized tract length: 36.1363636363636"
## [1] "Max lnL at grid value: 70"
## [1] "PSJS estimate: 66.9239527157328"

## [1] "In simulation, number of IGC events: 29 , mean potential tract length: 54.0344827586207 , mean realized tract length: 51.4137931034483"
## [1] "Max lnL at grid value: 25"
## [1] "PSJS estimate: 23.328324790887"

## [1] "In simulation, number of IGC events: 28 , mean potential tract length: 34 , mean realized tract length: 28.25"
## [1] "Max lnL at grid value: 45"
## [1] "PSJS estimate: 42.9131982178835"

## [1] "In simulation, number of IGC events: 23 , mean potential tract length: 61.304347826087 , mean realized tract length: 54.9130434782609"
## [1] "Max lnL at grid value: 25"
## [1] "PSJS estimate: 26.6500272671767"

## [1] "In simulation, number of IGC events: 33 , mean potential tract length: 60.5151515151515 , mean realized tract length: 52.2121212121212"
## [1] "Max lnL at grid value: 35"
## [1] "PSJS estimate: 32.7366328633833"

## [1] "In simulation, number of IGC events: 28 , mean potential tract length: 56.0714285714286 , mean realized tract length: 53.9642857142857"
## [1] "Max lnL at grid value: 45"
## [1] "PSJS estimate: 42.6751103103103"

## [1] "In simulation, number of IGC events: 28 , mean potential tract length: 41.3214285714286 , mean realized tract length: 39.5"
## [1] "Max lnL at grid value: 20"
## [1] "PSJS estimate: 17.5351242385874"

## [1] "In simulation, number of IGC events: 37 , mean potential tract length: 47.4054054054054 , mean realized tract length: 42.1621621621622"
## [1] "Max lnL at grid value: 40"
## [1] "PSJS estimate: 42.0679516107738"

## [1] "In simulation, number of IGC events: 29 , mean potential tract length: 47.8965517241379 , mean realized tract length: 47.551724137931"
## [1] "Max lnL at grid value: 25"
## [1] "PSJS estimate: 26.2594817529259"

## [1] "In simulation, number of IGC events: 16 , mean potential tract length: 43.5 , mean realized tract length: 36.5"
## [1] "Max lnL at grid value: 40"
## [1] "PSJS estimate: 41.4335015517569"

## [1] "In simulation, number of IGC events: 27 , mean potential tract length: 43.1851851851852 , mean realized tract length: 42.1111111111111"
## [1] "Max lnL at grid value: 15"
## [1] "PSJS estimate: 15.9165247984269"

## [1] "In simulation, number of IGC events: 27 , mean potential tract length: 51.6666666666667 , mean realized tract length: 51.6666666666667"
## [1] "Max lnL at grid value: 35"
## [1] "PSJS estimate: 34.9435587209407"

## [1] "In simulation, number of IGC events: 29 , mean potential tract length: 37.2758620689655 , mean realized tract length: 34.8620689655172"
## [1] "Max lnL at grid value: 30"
## [1] "PSJS estimate: 28.0988140272391"

## [1] "In simulation, number of IGC events: 31 , mean potential tract length: 70.6774193548387 , mean realized tract length: 64.7096774193548"
## [1] "Max lnL at grid value: 50"
## [1] "PSJS estimate: 48.4350491782479"

## [1] "In simulation, number of IGC events: 30 , mean potential tract length: 75.7 , mean realized tract length: 68.6666666666667"
## [1] "Max lnL at grid value: 70"
## [1] "PSJS estimate: 68.4325603274582"

## [1] "In simulation, number of IGC events: 31 , mean potential tract length: 56.7741935483871 , mean realized tract length: 43.3225806451613"
## [1] "Max lnL at grid value: 25"
## [1] "PSJS estimate: 24.9243491525757"

## [1] "In simulation, number of IGC events: 23 , mean potential tract length: 43.7826086956522 , mean realized tract length: 36.1304347826087"
## [1] "Max lnL at grid value: 25"
## [1] "PSJS estimate: 25.8560209795406"

## [1] "In simulation, number of IGC events: 30 , mean potential tract length: 44.1666666666667 , mean realized tract length: 39.3666666666667"
## [1] "Max lnL at grid value: 55"
## [1] "PSJS estimate: 55.8126823306154"

## [1] "In simulation, number of IGC events: 28 , mean potential tract length: 53.3571428571429 , mean realized tract length: 46.4642857142857"
## [1] "Max lnL at grid value: 40"
## [1] "PSJS estimate: 42.3980385298825"

## [1] "In simulation, number of IGC events: 32 , mean potential tract length: 67.25 , mean realized tract length: 56.78125"
## [1] "Max lnL at grid value: 20"
## [1] "PSJS estimate: 20.316523721997"

## [1] "In simulation, number of IGC events: 22 , mean potential tract length: 46.2272727272727 , mean realized tract length: 41.5454545454545"
## [1] "Max lnL at grid value: 35"
## [1] "PSJS estimate: 34.4787297004118"

## [1] "In simulation, number of IGC events: 19 , mean potential tract length: 55.1578947368421 , mean realized tract length: 52.5263157894737"
## [1] "Max lnL at grid value: 20"
## [1] "PSJS estimate: 17.6711339602354"

## [1] "In simulation, number of IGC events: 25 , mean potential tract length: 49.68 , mean realized tract length: 46.44"
## [1] "Max lnL at grid value: 30"
## [1] "PSJS estimate: 28.0551629559558"

## [1] "In simulation, number of IGC events: 29 , mean potential tract length: 47.3103448275862 , mean realized tract length: 47.3103448275862"
## [1] "Max lnL at grid value: 15"
## [1] "PSJS estimate: 14.2704366765976"

## [1] "In simulation, number of IGC events: 27 , mean potential tract length: 62.2962962962963 , mean realized tract length: 58.037037037037"
## [1] "Max lnL at grid value: 35"
## [1] "PSJS estimate: 36.4462673570927"

## [1] "In simulation, number of IGC events: 22 , mean potential tract length: 46.4545454545455 , mean realized tract length: 46.4545454545455"
## [1] "Max lnL at grid value: 25"
## [1] "PSJS estimate: 24.3071856409962"

## [1] "In simulation, number of IGC events: 29 , mean potential tract length: 61.448275862069 , mean realized tract length: 46.1034482758621"
## [1] "Max lnL at grid value: 30"
## [1] "PSJS estimate: 29.1833373597483"

## [1] "In simulation, number of IGC events: 24 , mean potential tract length: 41.0416666666667 , mean realized tract length: 38"
## [1] "Max lnL at grid value: 35"
## [1] "PSJS estimate: 37.0871546615615"

## [1] "In simulation, number of IGC events: 34 , mean potential tract length: 65.7352941176471 , mean realized tract length: 52.5882352941176"
## [1] "Max lnL at grid value: 35"
## [1] "PSJS estimate: 36.0507191972243"

## [1] "In simulation, number of IGC events: 20 , mean potential tract length: 51.8 , mean realized tract length: 47.45"
## [1] "Max lnL at grid value: 40"
## [1] "PSJS estimate: 39.7034674772249"

## [1] "In simulation, number of IGC events: 27 , mean potential tract length: 34.2592592592593 , mean realized tract length: 33.0740740740741"
## [1] "Max lnL at grid value: 70"
## [1] "PSJS estimate: 72.143177036402"

## [1] "In simulation, number of IGC events: 30 , mean potential tract length: 48.8666666666667 , mean realized tract length: 47.6666666666667"
## [1] "Max lnL at grid value: 55"
## [1] "PSJS estimate: 53.4838072854267"

## [1] "In simulation, number of IGC events: 26 , mean potential tract length: 44.7307692307692 , mean realized tract length: 40.3846153846154"
## [1] "Max lnL at grid value: 10"
## [1] "PSJS estimate: 10.8569661034689"

## [1] "In simulation, number of IGC events: 24 , mean potential tract length: 56.4166666666667 , mean realized tract length: 52.9166666666667"
## [1] "Max lnL at grid value: 40"
## [1] "PSJS estimate: 38.3706401635429"

## [1] "In simulation, number of IGC events: 22 , mean potential tract length: 44.4545454545455 , mean realized tract length: 43.5909090909091"
## [1] "Max lnL at grid value: 30"
## [1] "PSJS estimate: 32.2150543402945"

## [1] "In simulation, number of IGC events: 30 , mean potential tract length: 63.3333333333333 , mean realized tract length: 47.9666666666667"
## [1] "Max lnL at grid value: 60"
## [1] "PSJS estimate: 57.5435079308362"

## [1] "In simulation, number of IGC events: 29 , mean potential tract length: 60.2758620689655 , mean realized tract length: 55.9310344827586"
## [1] "Max lnL at grid value: 80"
## [1] "PSJS estimate: 76.7830049263824"

## [1] "In simulation, number of IGC events: 25 , mean potential tract length: 60.12 , mean realized tract length: 53.44"
## [1] "Max lnL at grid value: 25"
## [1] "PSJS estimate: 23.8283172833799"

## [1] "In simulation, number of IGC events: 22 , mean potential tract length: 42 , mean realized tract length: 34.4090909090909"
## [1] "Max lnL at grid value: 25"
## [1] "PSJS estimate: 22.5638082879767"

## [1] "In simulation, number of IGC events: 29 , mean potential tract length: 76.9655172413793 , mean realized tract length: 63.3103448275862"
## [1] "Max lnL at grid value: 45"
## [1] "PSJS estimate: 47.0249333210516"

## [1] "In simulation, number of IGC events: 24 , mean potential tract length: 49 , mean realized tract length: 42.75"
## [1] "Max lnL at grid value: 60"
## [1] "PSJS estimate: 61.6668889239587"

## [1] "In simulation, number of IGC events: 27 , mean potential tract length: 66.5555555555556 , mean realized tract length: 55.8518518518519"
## [1] "Max lnL at grid value: 60"
## [1] "PSJS estimate: 60.1502759117779"

## [1] "In simulation, number of IGC events: 32 , mean potential tract length: 36.75 , mean realized tract length: 36.3125"
## [1] "Max lnL at grid value: 30"
## [1] "PSJS estimate: 31.6736418888907"

## [1] "In simulation, number of IGC events: 28 , mean potential tract length: 40.3928571428571 , mean realized tract length: 37.8214285714286"
## [1] "Max lnL at grid value: 25"
## [1] "PSJS estimate: 25.4875109510536"

## [1] "In simulation, number of IGC events: 29 , mean potential tract length: 44.9655172413793 , mean realized tract length: 39.3793103448276"
## [1] "Max lnL at grid value: 35"
## [1] "PSJS estimate: 34.2500258124985"

## [1] "In simulation, number of IGC events: 30 , mean potential tract length: 38.0666666666667 , mean realized tract length: 33.6666666666667"
## [1] "Max lnL at grid value: 25"
## [1] "PSJS estimate: 24.6238690643698"

## [1] "In simulation, number of IGC events: 32 , mean potential tract length: 45.28125 , mean realized tract length: 42.75"
## [1] "Max lnL at grid value: 50"
## [1] "PSJS estimate: 49.5684283066384"

## [1] "In simulation, number of IGC events: 29 , mean potential tract length: 50.9310344827586 , mean realized tract length: 38.3793103448276"
## [1] "Max lnL at grid value: 100"
## [1] "PSJS estimate: 105.797229994723"

## [1] "In simulation, number of IGC events: 27 , mean potential tract length: 64.0740740740741 , mean realized tract length: 45.2592592592593"
## [1] "Max lnL at grid value: 35"
## [1] "PSJS estimate: 35.6807598166827"

## [1] "In simulation, number of IGC events: 26 , mean potential tract length: 44.1153846153846 , mean realized tract length: 42.5"
## [1] "Max lnL at grid value: 35"
## [1] "PSJS estimate: 36.5954718933644"

## [1] "In simulation, number of IGC events: 28 , mean potential tract length: 53.4642857142857 , mean realized tract length: 36.5714285714286"
## [1] "Max lnL at grid value: 25"
## [1] "PSJS estimate: 26.2638231818596"

## [1] "In simulation, number of IGC events: 24 , mean potential tract length: 47.0833333333333 , mean realized tract length: 41.875"
## [1] "Max lnL at grid value: 30"
## [1] "PSJS estimate: 30.1654131811983"

## [1] "In simulation, number of IGC events: 29 , mean potential tract length: 49.8620689655172 , mean realized tract length: 45.2758620689655"
## [1] "Max lnL at grid value: 60"
## [1] "PSJS estimate: 59.0176461808477"

## [1] "In simulation, number of IGC events: 25 , mean potential tract length: 47.36 , mean realized tract length: 44.6"
## [1] "Max lnL at grid value: 25"
## [1] "PSJS estimate: 23.7181952490129"

## [1] "In simulation, number of IGC events: 25 , mean potential tract length: 52.28 , mean realized tract length: 50.08"
## [1] "Max lnL at grid value: 60"
## [1] "PSJS estimate: 64.0596010859126"

Save workspace

save.image("./GridWeight.RData")